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1 Introduction and definitions 



1.1 Introduction 

Semi-infinite lattice systems with free surface were studied extensively by 
mean field and numerical approaches. They exhibit a rich variety of critical 
behaviour depending on the relative strengths of the surface and bulk cou- 
pling constants. Ordinary critical behaviour occurs when the surface spins 
order at the bulk critical temperature as the temperature is lowered. When 
the surface couplings are sufficiently enhanced, the surface spins order above 
the bulk critical temperature in a transition known as the surface transition. 
As the temperature is then lowered through the bulk critical temperature 
the extraordinary transition takes place. Finally, a special critical behaviour 
arises if the couplings are adjust such that the bulk and surface transitions 
coincide. We refer the reader to Binder jHj for a review on the subject. Exact 
results may be found in pp, [2E|- A rigorous approach of the semi- infinite 
Ising model is proposed in a series of articles by Frohlich and Pfister, [[•*)]- 
[T7j . In particular, they observed that the surface transition does not occur 
for temperatures less than the critical bulk temperature (see Figure 1 in |15j). 

A new scheme, was found by Lipowsky in the case of the g-state Potts 
model. In the many component limit q — * oo, both mean field theory and 
Migdal-Kadanoff renormalization group predict a new phase where the bulk 
is ordered while the free surface is disordered. The relative strengths of the 
surface and bulk couplings can be chosen such that when the temperature 
is lowered, first the bulk spins order through an extraordinary transition 
and then surface spins order through a surface transition, or a new special 
transition arises (see diagram Figure 1 in [26]). Intuitively, one can think that 
the low temperature ordered bulk acts upon the spin on the boundary surface 
like an external magnetic field. Roughly speaking, one can imagine that this 
effect modifies the Hamiltonian in the boundary surface by a term containing 
a magnetic field. Then for the Ising case a surface phase transition cannot 
occur when the bulk is ordered. However, this observation is not true for 
the q-state Potts model where it is known for q large that a disordered state 
arises even in a presence of a magnetic field and coexists with the ordered 
one aligned with the magnetic field |5J. 

The aim of this work is to prove that for q large the mean field or renor- 
malization group prediction are indeed correct . We will study both surface 
phase transitions (of the semi-infinite g-state Potts model for bulk dimension 
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d > 3), using the Fortuin-Kasteleyn representation of this model. In doing 
so we think that the method can be extended to a broader class of systems. 

Let us sketch out the main scheme of the approach. Since we have to 
study the surface phase transitions and thus the possible coexistence of "sur- 
face phases" (that actually holds), we define following [T3] two surface free 
energies corresponding respectively to the high and low temperature bulk. 
Here by "surface phases" we mean states invariant under horizontal trans- 
lations. The ratio of partition functions defining the surface free energies 
is then expanded using high or low convergent expansion. This leads to 
horizontally translation invariant systems called hydra models that are then 
analyzed with the help of Pirogov-Sinai theory (22 . Actually, the theory de- 
veloped for the study of translation invariant bulk states can be implemented 
to the case under consideration. 

The paper is organized as follow. The surface free energies are defined in 
the next subsection. In Section 2, we present the Fortuin-Kasteleyn represen- 
tation and the ground state analysis of the model. Then we give convergent 
cluster expansions of the bulk partition function in the high and low (bulk) 
temperature regime. The hydra system for the high temperature bulk case 
is presented in Subsection 3.1. Subsections 3.2 and 3.3 are devoted to the 
large q analysis of the hydra model. Finally, some proofs are postponed to 
the appendix. 

1.2 Surface free energies 

Consider a ferromagnetic Potts model on the semi-infinite lattice L = Z d_1 x 
Z + of dimension d > 3. At each site i = {ii,---,id} £ L, with i a 6 Z for 
a = 1, d— 1 and id G Z + , there is a spin variable <7j taking its values in the 
set Q = {0, 1, . . . , q — l}. We let d(i,j) = max a=lr . ^ d \i a — j a \ be the distance 
between two sites, d(i,Q) = mmj^n d(i, j) be the distance between the site 
i and a subset Q C L, and d(Q,Q') = minj g n,jeo' d(i, j) be the distance 
between two subsets of L . The Hamiltonian of the system is given by 



where the sum runs over nearest neighbor pairs (i.e. at Euclidean 

distance d^(i,j) = 1) of a finite subset Q C L, and S is the Kronecker 
symbol: 5(<Ji,aj) = 1 if G{ = <jj, and otherwise. The coupling constants 




(1.1) 



3 



Kij can take two values according both i and j belong to the boundary layer 
L = {i G L | id = 0}, or not: 

K«={ K , "a CL ° (1.2) 

I J > U otherwise 

The partition function is defined by: 

^(fi) = 5>-^Vn (1-3) 

Here the sum is over configurations cr^ G Q Q , (3 is the inverse temperature, 
and Xn is a characteristic function giving the boundary conditions. In par- 
ticular, we will consider the following boundary conditions: 

• the free boundary condition: xh — 1 

• the ordered boundary condition: Xn — Yli^an &{ a i-> 0); where the bound- 
ary of Q is the set of sites of Q at distance one to its complement 
Oft = {i G Q : d(i,L\fi) = 1}. 

• the free boundary condition in the bulk and ordered boundary condition 
on the surface: Xn — Ylied n &{ a ii 0)' where d s Q = dQ n L 

• the ordered boundary condition in the bulk and free boundary condition 
on the surface: Xn = Ylied b n 0),where d^l = dfl \ d s Q. 

Let us now consider the finite box 

Q = {i G L | max \i a \ < L, < i d < M} 

a=l,...,d— 1 

its projection E = Q D L = {i G Q \ id = 0} on the boundary layer and its 
bulk part A = = {i G tt | 1 < i d < M}. 

Following Frohlich and Pfister j5 3 we introduce two surface free energies, 
the free (or disordered) surface free energy with free boundary condition in 
the bulk, 

1 Z i (Vt) 
Qi = — lim — — - lim In——— (1-4) 

y L^oo |E| M-,00 Qf(A) V ^ 

and the ordered surface free energy corresponding to ordered boundary con- 
dition in the bulk, 

1 , , z°(n) . , 

g Q = — lim 7=rr lim In——— (1.5) 

^ L^oo |E| M^oo Q°(A) V ' 



Here |E| = (2L + l) d 1 is the number of lattice site in £, and Q f (A) and 
Q°(A) are the following bulk partition functions: 

Q f (A) = ]Texp{/3J S (°i>°i)} (1-6) 

(*,i)cA 

Q°(A) = ^exp{/3J £ J] %*.°) ( L7 ) 

(i,j)cA i&dA 

where the two sums are over configurations a\ G Q A . The surface free ener- 
gies do not depend on the boundary condition on the surface, in particular 
one can replace Z { {Vt) by Z fo (fi) in JO} and replace Z°(Q) by Z o{ {Vt) in 
()1.5p . The partial derivative of the surface free energy with respect to j3K 
represents the mean surface energy. As a result of this paper we get, for q 
large and e^ J — 1 < q 1 ^, that the mean surface energy g^gi is discontinuous 

near /3K = 1n (l + g 1 /^ 1 )). 

Namely, let (-) p denote the infinite volume expectation corresponding to 
the boundary condition p: 

(f) p (/3J,f3K) = lim V fe-^xl 

defined for local observable / and let e~ T be defined by (|3.22|) below. As 
a consequence of our main result (Theorem 13.61 in Section 3), we have the 
following 

Corollary 1.1 Assume that e@ J — 1 < q x l d and q is large enough, then there 
exists a unique value K t such that 

(i) for any n.n. pair ij of the surface 

(SfaajWiPJrfK) <0(e~ T ) for K <K t 
(5(a 2 ,,a J ))^J,/?if)>l-0(e- T ) for K > K t 

(ii) for any n.n. pair ij between the surface and the first layer 

{5{a u a j )Y{pj^K)<0{e- T ) 
(S(<T il <T j )y°(pj,pK) <0(e~ T ) 
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In that theorem the ratios of the partition functions entering in the defini- 
tion of the surface free energy (with both Z { (Q) and Z fo (Q)) are expressed 
in terms of partition functions of gas of polymers interacting through a two- 
body hard-core exclusion potential. For e@ J — 1 < q l / d and q large, the 
associated activities are small according the values of K namely for K < K t 
with the free boundary condition and for K > K t with the free-ordered 
boundary condition. The system is then controlled by convergent cluster 
expansion. 

In the next subsection, we give the expression of the Fortuin-Kasteleyn 
representation of the partition function Z P (Q). We also present, in Theo- 
rems 2.1 and 2.2, the high and low temperature expansions of the Fortuin- 
Kasteleyn representation of the bulk partition function Q(A) that converge 
respectively for e^ J — 1 < q 1 ' d and e^ J — 1 > q l ^ d whenever q is large enough. 

2 Random cluster models 

2.1 The Fortuin-Kasteleyn (FK) representation 

By using the expansion e /3KijS ^ i ' a ^ = 1 + {eP Kli — l)5((7j, <7j), we obtain the 
Fortuin-Kasteleyn representation ^3] of the partition function: 

Z?(Q) = II ^ Kl3 - l ^ K(X) ( 21 ) 

xcB(n) (i,j)ex 

where B(Q) = : i G f2, j G Q} is the set of bonds with both endpoints 

belonging to Q, and N^(X) is the number of connected components (regard- 
ing an isolated site i G Q as a component) of a given X C B(Q). These 
numbers depend on the considered boundary condition; introducing S(X) as 
the set of sites that belong to some bond of X and C(X\V) as the number 
of connected components (single sites are not included) of X that do not 
intersect the set of sites V, they are given by: 



<P0 = 


|0| 


-\S(X)\ + C(X\V>) 




\n\ 


- \S(X) Udn\ +C(X\dQ) 


nL°(x) = 


\n\ 


- \S(X) Ud s Q\ +C{X\d s Q) 




\n\ 


- \S(X)Ud b Q\+C{X\d b fl) 
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Hereafter \E\ denotes the number of elements of the set E. The Boltzmann 
weight in (|2.1|) can be written in terms of the following Hamiltonian: 

H*(X) = -P.\X.\ - (3 b \X b \ - N*(X) + p\ (2.2) 

where X s = X n B(L ), X b = X \ X s , and 

ln(e^ - 1) 



mq 

_ ln(e^ J - 

Pb 



(2.3) 



In q 

Actually, one has: 

2?(n) = qM J2 ( 2 - 4 ) 

XcB(Q) 

Hence from this representation where the configurations are given by sub- 
sets of bonds of the lattice, we are left to analyze a model with an inverse 
temperature given by \nq and Hamiltonian given by (j2.2j) where (3 S and j3 b 
are independent (real-valued) coupling constants. The next subsection is 
devoted to the study of the ground states of this model. 



2.2 Ground states 

To study the diagram of ground states of the Fortuin-Kasteleyn representa- 
tion, it is convenient to introduce the formal Hamiltonian 

H(X) = -f3 s \X s \ - (3 b \X b \ + \S(X)\ - C(X) (2.5) 

where C(X) = C(X|0) is the number of connected components of X, and 
the relative Hamiltonian 

H(Y\X) = H(Y) - H(X) (2.6) 

Note, that this last expression makes sense for any configurations X and Y 
that coincide almost everywhere (a.e.), i.e. that differ only on a finite set 
of bonds. A configuration X is called ground state, if for any configuration 
Y = X a.e., one has 

H(Y\X) > 

For a fixed bond b of the lattice let n b be the occupation number: n b — 1 
if b belongs to X and otherwise. The structure of ground states invariant 
under horizontal translations is indicated in Fig. 1 (for positive couplings j3 s 
and f3 b ). 
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FIGURE 1 : The diagram of ground states. 
There are four regions with a unique ground state: 

1. In region I ( < i , (3 S < -^j, the ground state is the configuration 
.V . 

2. In region II ((3 b < | , (3 S > it is the configuration X fo = B(h ). 

3. In region III ( f3 h > | , f3 s < ^-(1 — /3b), it is the configuration X of = 
5(L\L ). 

4. In region IV ( > | , (3 S > ^zj(l — /9b), it is the configuration X° = 
£(L). 

On the separating lines between any two regions the corresponding ground 
states coexists. The g.s. X { coexists with X {o and X° at the point Si and 
with X of and X° at the point S^- 

On the part of the line fid = 1/d above and on Si, there are furthermore 
infinitely many g.s. that coexist with the previous ones, namely all configu- 
rations say Xk, k = 2,3, with occupied bonds rib — 1 below and on the 
k layer and empty bonds = otherwise. Analogously, on the part of the 
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line (3d = 1/d below and on S2 there are also infinitely many g.s. coexisting 
with the previous ones, namely all configurations X' k , k = 2,3, with oc- 
cupied bonds above and on the k layer and empty bonds otherwise. Finally, 
when f3 s = and /3& > 1, comes in addition to this list the configuration 
X oi = B(L) \B(Lo). 

This suggest that for large q the extraordinary transitions will occur 
through a sequence of layering transitions. Here we concentrate on the sur- 
face transitions, but let us mention a particular study of layering transitions 
and entropic repulsion for a simplified Solid-On-Solid type model in 0, [7j. 

2.3 FK representation in the high temperature regime 

In terms of the Fortuin-Kasteleyn representation, the partition function 
<2 f (A) reads: 

Q\A)=q^ q MYHS{¥)l+C{Y) (2.7) 

YcB(A) 

We will first recall that this system is described by a gas of polymers inter- 
acting through two-body hard-core potential and with activity given by (J2.8|) 
below. 

Indeed, first let us call polymer, a connected subset of B(h) (in the M, d 
sense). We will use V(A) to denote the set of polymers whose bonds belong 
to -B(A). Two polymers 71 and 72 are said compatible (we will write 71 ~ 72) 
if they do not intersect, i.e. if the set 5(71) R £(72) is empty. They are said 
incompatible otherwise and we will use the notation 71 ^ 72. A family of 
polymers is said compatible if any two polymers of the family are compatible 
and we will use P(A) to denote the set of compatible families of polymers 

7 e P(A). 

Clearly each subset Y of -B(A) may be decomposed (in a unique way) in 
a compatible family of polymers Y = {71, . . . ,7 n }- Introducing the activity 
of polymers by 

^=^171-15(7)1+1 (2. 8 ) 

one has: 

Q f (A)=g |A| Yl Il^fr) ( 2 - 9 ) 

yeP(A) 7sv 

(with the sum running over compatible families of polymers including the 
empty-set with weight equal to 1). 
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We will now introduce multi-indexes in order to write the logarithm of this 
partition function as a sum over these multi-indexes (see |27j). A multi-index 
C is a function from the set V(A) into the set of non negative integers, and 
we let suppC = {76 'P(A) : C(-y) > 1}. We define the truncated functional 

where the factor a(C) is a combinatoric factor defined in terms of the con- 
nectivity properties of the graph G(C) with vertices corresponding to 7 G 
supp C (there are C (7) vertices for each 7 G supp C) that are connected by 
an edge whenever the corresponding polymers are incompatible). Namely, 
a(C) = and hence $o(C) = unless G{C) is a connected graph in which 
case C is called a cluster and 

a(C)= (-l) |e(G)l (2.H) 

GCG{C) 

Here the sum goes over connected subgraphs G whose vertices coincide with 
the vertices of G(C) and |e(G)| is the number of edges of the graph G. If the 
cluster C contains only one polymer, then 0(7) = 1. In other words, the set 
of all bonds of polymers belonging to a cluster C is connected. The support 
of a cluster is thus a polymer and it is then convenient to define the following 
new truncated functional 

$( 7 )= Yl ( 2 - 12 ) 

C:supp C=-f 

As a result we have 

Theorem 2.1 Assume that fa < 1/d and c^v^q'^^ < I, where v& = (2d) 2 
and c = (a/3 + 2)e^+T ; then 

Qf(A) = g |A| exp J £ $(7) 1 (2.13) 

with a sum running over (non-empty) polymers of A, and the truncated func- 
tional $ satisfies the estimates 

< ItI fco^g-^V 7 ' ( 2 - 14 ) 
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Proof. We first notice that easy geometrical observations [231 121] gives 
the inequality ^ — |5(7)| + 1 < which in turn implies: 

Hi) < q-G-MW (2.15) 

On the other hand the number of polymers 7 of length I7I = n containing a 
given site is less than (vd) n - This follows from the fact that for any polymer 
there exists a path going through the center of every bond at most twice (see 
jSUJ). Then we have X) 73a; <Pf(j) < 00 as soon as v d q~d + P b < 1. 



According ref. [2I|, under the condition 



Vf(7o) < (e M(7o) - 1) 



where \i is a positive function, one has 



exp 



5>(7) 



7^70 



(2.16) 



veP(A)7ev c 

where the last sum runs over (non empty) clusters of polymers in V(A), and 
the truncated functional $ satisfies the estimate 

E l $ o(^)|<M7i) (2.18) 

C:C( 7 i)>l 

We choose ^(7) = (az^)~' 7 ' to get by taking into account the above remark 
on the entropy of polymers 

00 n~ l 2 

E A*(7) < |5(7o)| E^" ^ (iTol + 1)t— :=i < 73 7^1 (2-19) 

79^70 71=1 

for a > 1. Since /i(7o) < e^ 70 - 1 — 1, the bound (|2.15|) gives that the condition 
(|2.16|) will be satisfied whenever 

u d q-^ +l3b ae^ < 1 (2.20) 

The value ao = \/3 + 2 that minimizes the function ae^ 31 provides the con- 
ditions given in the theorem. By definition (|2.12|1 . the relation ()2.17j) yields 
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(jHHD while the relations (|2~TH|) and (gUD yield: 

|$(7o)|<e" A ^l ]T |$o(C)ine AC(7) (2-21) 

C:supp C=7o 7 

< e " A|7 °'E E l*o(C)ine AC,(7) <^-rl7o|e- A ^l 



7^70 C:C(7)>1 



a - 1 



provided z/^e A g~^ + ^ i, aoe"o- T < 1; we choose for A the value that realizes the 
equality to get (by taking into account that ^-j- < 1) the estimate (j2.14|) . ■ 
Let us mention that a recent study of these expansions for the Potts 
model on a general graph, i.e. not restricted to nearest neighbor interactions, 
is given in ref. 



2.4 FK representation in the low temperature regime 

We now derive an expansion of the bulk partition function Q°(A) at "tem- 
perature" (3b > \- The FK representation of this partition function reads 

Q°(A) = ^ q MYl+N * {Y) = q MB(A)l ^ q~ MB ^ Yl+N ^ Y) (2.22) 

KcB(A) YcB(A) 

where N£(Y) = \A\-\S(X) U <9A| + C(X|<9A). The expansion is mainly based 
on a duality property and we first recall geometrical results on Poincare and 
Alexander duality (see e.g. |25|.|2|.P1.|I3|). 

We first consider the lattice Z d and the associated cell-complex L whose 
objects s p are called p-cells (0 < p < d): 0-cells are vertices, 1-cells are bonds, 
2-cells are plaquettes etc.: a p-cell may be represented as (x; <T\e\, a p e p ) 
where x G Z d , (ei, e<i) is an orthonormal base of M. d and a a = ±l,a = 
1, d. Consider also the dual lattice 

r i i 

(Z d )* = < x = (xi + -, x d + -) : x a e Z, a = 1, d 

and the associated cell complex L*. There is a one to-one correspondence 

s P ~ s* d _ p (2.23) 

between p-cells of the complex L and the d — p-cells of L*. In particular to 
each bond s\ corresponds the hypercube s^-i that crosses si in its middle. 
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The dual E* of a subset E C L is the subset of elements of L* that are in 
the one-to-one correspondence (|2.23|) with the elements of E. 

We now turn to the Alexander duality in the particular case under con- 
sideration in this paper. Let Y C B(A) be a set of bonds. We define the 
A-dual of Y as 

Y = (B(A)\Y)* (2.24) 
As a property of Alexander duality one has 



Y 



|B(A)\y| (2.25) 



Nl(Y) = N d (Y) (2.26) 

where N C \(Y) denotes the number of independent closed (d — l)-surfaces of 
Y. We thus get 

Q°(A) = q^ B ^\ q~ M?l+Nd{9) (2.27) 

rc[-B(A)]* 

Now as in the previous subsection, the system can be described by a gas 
of polymers interacting through hard core exclusion potential. Indeed, we 
introduce (hat)-polymers as connected subsets (in the M. d sense) of (d — 1)- 
cells of L* and let V(A) denote the set of polymers whose (d— l)-cells belong 
to [5(A)]*. Two polymers 71 and 72 are compatible (we will write 71 ~ 72) if 
they do not intersect and incompatible otherwise (we will write 71 00 72). A 
family of polymers is said compatible if any two polymers of the family are 
compatible and we will use P(A) to denote the set of compatible families of 
polymers 7 G V(A). Introducing the activity of polymers by 

(p tf) = q -Mi\+N al (j) (2>28) 

one has: 

Q o {A)=q MB(A)\ ^ JJ^) (2 29) 

YeP(A)jeY 

(with the sum running over compatible families of polymers including the 
empty-set with weight equal to 1). As in the previous subsection, we introduce 
multi-indexes C as functions from the set ^(A) into the set of non negative 
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integers, and the truncated functional 



a{C) 



]>°(7) 



cm 



8( 7 ) = Yl 



C:supp C=7 



(2.30) 
(2.31) 



Here suppC = E V(A) : (^(7) > l| and a(C) is defined by (l2~TTJ) with a 

graph G(C) with vertices corresponding to 7 G supp C that are connected by 
an edge when the corresponding polymers are incompatible. Again, $o(C) 
( ) if ( ' is not a cluster i.e. if the set of all d — 1-cells of polymers belonging 
to supp C is not connected. 



Theorem 2.2 Assume that (3^ > 1/d and Co^g 



d 2 2 4(d-i) ! and c = 1 + 2 d " 2 (l + VT+2F 1 ) 



cxp 



< 1, where % 



, then 



Q°{A) = e 



0b\B(A)\ 



76P(A) 



(2.32) 



with a sum running over (non-empty) polymers, and the truncated functional 
$ satisfies the estimates 



$(7) 



< 7 Wdq 



l\l7l 



^2.33) 



Proof. We first notice that the obvious geometrical inequality N c ifi) < 
\j\/d implies: 



^0(7) < q-^-^ 1 



[2.34) 



On the other hand the number of polymers 7 of length I7I = n containing a 
given vertex is less than (pd) n - To see it, we first observe that the number of 
p-cells that share a same vertex equals 2 P (^) : 2 P is the choice for the signs of 

the ex's and the binomial coefficient (^) is the choice for p vectors e a among 
(ei, ed). This implies (by duality) that a (d— l)-cell contains 2 d ~ 1 vertices 
and hence that the number of (d — l)-cells connected (i.e. sharing at least a 
vertex) with a given {d — l)-cell is less than 2^ _1 ^(d2^ -1 ^ — 1). Finally, one 
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uses the fact that for any polymer of length n there exists a path of length 
less than 2n going through the center of every (d — l)-cell. 
For a function = (aP^) - '^' we get 

^ ^ -1 d L 

77^70 71=1 

if a > 1 and the convergence condition (j2.15|) for the case under consideration 
will be satisfied whenever 

9 d q' Ph+ ^ae^ < 1 (2.36) 

Here we choose a = 1 + 2 d ~ 2 (l + y/l + 2 3 ~ d ) that minimizes the function 

ae a ~ 1 . 

The remainder of the proof is analog to that of Theorem 12.11 ■ 

3 Surface transition in the high temperature 
bulk regime 

3.1 Hydra model 

To study the surface free energy g>f we shall express the ratio of partition 
functions entering in its definition (jl.4j) in terms of a partition function of 
geometrical objects to be called hydras. 

Definition 3.1 A connected set of bonds 5 C B(Q) is called hydra if there 
exists a bond of 5 which has one of its endpoint on the boundary surface S 
(i.e. ifS(5) nS ^ 0;. 

Any bond of the hydra with one endpoint in the boundary surface S and 
one endpoint in the bulk A is called neck of the hydra. 

A connected component of bonds of the hydra with two endpoints in the 
boundary surface is called foot of the hydra. 

A connected component of bonds of the hydra with two endpoints in the 
bulk is called head of the hydra (heads of hydras are polymers introduced in 
Subsection 2.3). 
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FIGURE 2 : A hydra, in two dimensions (a dimension not considered in this paper), 
with 2 feet (components of full lines), 5 necks (dotted lines), and 3 heads (components of 

dashed lines). 

We let 7i(Q) denote the set of hydras of Q. Two hydras 5\ and 62 are said 
compatible (we will write 5\ ~ 82) if S(5i) H Sfa) 7^ 0. A family of hydras 
is said compatible if any two hydras of the family are compatible and we let 
H(f2) denote the set of compatible families of hydras 5 G TC(fl). 

Clearly, a connected subset of bonds included in B(Q) is either a hydra 
5 G H(Q) or a polymer 7 G P(A) (defined in subsection 2.3) and thus any 
subset of B G B(Q) is a disjoint union of a compatible family of hydras 
X G H(fi) with a compatible family of polymers Y G P(A). 

Then, the partition function Z P (Q) defined by (j2.4j) reads 

z^) = q w r mx) E 11^) 

JfeH(Sl) YeP{A):Y~X -y&Y 

where we take into account that the Boltzmann weight q~ H n( Y ) of the subsets 

Y C B(A) coincides with the bulk one given in ()2.7|) and the compatibility 

Y ~ X means S(Y) n S(X) = 0. 

According to Subsection 2.3, the last sum in the RHS of the above formula 

can be exponentiated as: exp |2^ 7 eP(A)- 7 ~x ^(t)|- Hence dividing the above 
partition function by the partition function Q f (A) we get by taking into 
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account Theorem 12.11 



s"(n) 



Z p {VL) 



(3.2) 



= q 



xeH(^) 




7£P(A);7^X 



Hereafter the incompatibility 7 X means 5(7) fl S'(X) 7^ 0. 

H P (Q) is thus the partition function of a gas of hydras X = {81, . . . , 5 n } 
whose activities are defined from the Hamiltonian Hf l {X); they interact 
through hard-core exclusion potential and through a long range interaction 
potential (decaying exponentially in the distance under the hypothesis of 
Theorem 12. 1J) defined on the polymers of the bulk. 

Notice that if we neglect this long range potential, the system of hydras 
will reduce itself to a (d — 1) Potts model (in the FK representation), when 
we moreover restrict to consider only hydras without neck and without head. 
As well known, the (d — 1) Potts model undergoes a temperature driven 
first order phase transition (at some inverse temperature ~ (d — l)~ 1 lnq), 
whenever q is large enough and d > 3. We will show that it is also the 
case for the hydra model ()3.2j) implementing the fact that the necks and the 
heads of hydras modify only weakly their activities and that the long range 
interaction potential decays exponentially (the needed assumptions are close 
to those of Theorem 12. 1|) . 

To this end it is convenient to first rewrite this potential in terms of a 
model of aggregates, a technique already developed in many analog situations 
(see e.g. [TS].[TTj). Let us introduce the (real- valued) functional 



defined on polymers 7 G V(A). An aggregate A is a family of polymers 
whose support, supp A = U 7£j 47, is connected. Two aggregates A\ and A 2 
are said compatible if supp A\ fl supp A 2 = 0. A family of aggregates is said 
compatible if any two aggregates of the family are compatible and we will 
use A (A) to denote the set of compatible families of aggregates. Introducing 
the statistical weight of aggregates by 



*(t) 



= e 



*(7) 



1 



(3.3) 



(A) = n *( 7 ) 



(3.4) 



7GA 
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we then get: 



exp 



e n 



7S"P(A); 76^(A); 



e n^) ( 3 - 5 ) 



i£A(A) 

where A ^ X means that every polymer of the aggregate A is incompatible 
with X. Since the support of aggregates is a connected set of bonds i.e. a 
polymer, it is convenient (as it was done for clusters in Subsection 2.3) to 
sum the statistical weights (|3.4|) over aggregates with same support. We thus 
define the weight 

A:supp A=7 

with A oo X, to get 

2 p (fi)=g |S| E ^ W E 11^(7) (3-7) 

xeH(Q) yeP(A) -y& 

The system is thus described by two families: a compatible family of 
hydras and a compatible family of polymers each of these polymers being 
incompatible with the family of hydras. 

The next subsections are devoted to the large q analysis of this model in 
the frame of Pirogov-Sinai theory. As well known in this theory, it is useful 
to introduce diluted and crystal partition functions for which an important 
set of recurrence relation holds (see Lemma 3.2 below). Let us introduce the 
diluted partition functions for our system . 

We define the envelope E(M) of a set of sites of the lattice L as the set 
of bonds of B(h) that contain an endpoint in M. Next it is convenient to 
define local observables with which the main terms of the Hamiltonian (the 
number of bonds and the number of sites) can be written. 

Let the energy per site of a configuration X C B(h) be: 

e t (X) = \E(i) nx s |-^ \E{i) n X b \ + \i n S{X)\ (3.8) 
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where X s = XnB(L ) and X b = X\X S . We define, for (any) volume ficL, 
the diluted Hamiltonian of a configuration X = X p 

H$\X) = d( X ) ~ C{X) + C{X?) (3.9) 

where the number of connected components C(X f ) of the ground state X f 
equals 0, and the number of connected component C(X fo ) of the ground state 
X io equals 1. 

Note that even though our model is defined on a <i-dimensional box Q it 
has a (d — l)-dimensional structure and the highest order of the logarithm 
of partition functions behaves like 0(|E|). It will be convenient to consider 
VI as a set of lines. We let a line L(x) be a cylinder set of sites of L whose 
projection on the boundary layer is the site x and whose height is less than a 
given number M: L(x) = {i E L ...,id-i) = x G L ,i rf < M}. We let h M 
be the set of all such lines. For Q C I^m, we let E = f2nL , be its projection 
on the boundary layer, A = Q \ E and = |E| be the number of lines 
of f2. We will use W(il) to denote the set of compatible families of hydras 
defined on the envelope E(Q) that coincide with the ground state X p on the 
envelope E(dQ) and on the envelope E(Lm \ Q), and use P dl1 (A) to denote 
the compatible families of polymers defined on E(Q) \ (E(E) U E(dQ)). 

The diluted partition function is defined by 

-f(^)= e ^ (x) e n^) ( 3 - iQ ) 

xeHP(Ji) y e p dil (A) T eY 

Notice that the diluted Hamiltonian on ground states reads on set of lines 

Q C L M : 

H™(X*)=e p \\n\\ (3.11) 

where 

e fo = i-(d-m (3 - 12) 

Up to a boundary term 0(<9E) one has lnS p (Q) = lng + lnS^ 1 ' (fl) so 
that 

-hm^LlnSf (n)= fl t + ln (Z (3.13) 

where Q | L means that we take first the limit M — > oo and then the limit 
E | L in the van-Hove or Fisher sense. 
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Notice also the following bound 



°-^: lnS " il(fi) - l|fi|l(rf_1)lng (3 - 14) 



3.2 contours 

Let Q C Lm and (X, Y) be a configuration of our system in Q: X G 
H. p ({l),Y G P dil (fi),F o« X. A site i G fl is called p-correct, if X coin- 
cides with the ground state X p on the bonds of the envelope E(i) and Y 
does not contain any bond of E(i). A line is called p-correct if all the sites 
of the line are p-correct. In other words a line L is f-correct if X n E(L) = 
and Y n = 0; a line L is fo-correct if X n = 5(L ) n £(L) and 

Y H = 0. Sites and lines that are not p-correct are called incorrect. 

The set of incorrect lines of a configuration (X, F) is called boundary of 
the configuration (X,Y). 

A triplet T = {supp T, X(r), F(r)}, where supp T is a maximal connected 
subset of the boundary of the configuration (X, Y) called support of T, X(T) 
the restriction of X to the envelope -E'(suppr) of the support of T, and 
Y(T) the restriction of Y to i?(supp T), is called contour of the configuration 
(X, Y). Hereafter a set of sites is called connected if the graph that joins all 
the sites i,j of this set with d(i,j) < 1 is connected. 

A triplet T = {supp T, X(r), F(r)}, where suppT is a connected set of 
lines is called contour if there exists a configuration (X, Y) such that T is a 
contour of (X,Y). We will use |T| to denote the number of incorrect points 
of supp T and ||r|| to denote the number of lines of supp T . 

Consider the configuration having T as unique contour; it will be denoted 
(X r , Y r ). Let L p (r) be the set of p-correct lines of \ suppT. Obviously, 
either a component of Lf(T) is infinite or a component of Lf (T) is infinite. 
In the first case T is called contour of the free class or f-contour and in 
the second case it is called fo-contour. When T is a p-contour (we will 
let T p denote such contours) we use Ext T to denote the unique infinite 
component of L p (T); this component is called exterior of the contour. The 
set of remaining components of L p (T) is denoted Int p r and the set L m ^ p {T) 
is denoted Int m r. The union Intr = IntfL U Intf o r is called interior of the 
contour and V(r) = supp T U Intr. 

Two contours r\ and T 2 are said compatible if the union of their supports 
is not connected. They are mutually compatible external contours if V(ri) C 
Ext T 2 and V(T 2 ) C Ext r\. 
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We will use G(T P ) to denote the set of configurations having T p as unique 
external contour. The crystal partition function is then defined by : 

E CT (T P ) = J2 q' H ^ iX) JJVCt) (3-15) 

(X,Y)eG(TP) -y&Y 

Lemma 3.2 The following set of recurrence equations holds : 

n 

~f(fi)= Yl q- epllExtll U ECr ^ ( 3 - 16 ) 

{rf,...,rR}«* i=i 

i/ere i/ie sum is over families {T p , . . . , rP} e:ri of mutually compatible external 
contours in Q (suppY p C f2j„t = {i G fi : d(i,L^ \ O) > ||-Krf|| = ||fi|| — 
^2\\V(T P )\\ where ||^(rf)|| is the number of lines ofV(T p ); 

i 

S cr (P0 = g{T p ) 11 Z%(Int m T p ) (3.17) 

me{fJo} 

where: 

Proof. We have only to observe that for any X e H p (f2) 

AnW = ^ i^prPn + ^e p ||L P (X) n Q\\ (3.19) 
r P 

where the sum is over all contours of the boundary of the configuration 
(X, Y = 0) and ||L p (X) PI is the number of p-correct lines inside Q of this 
configuration. ■ 

Lemma f3.2l gives the following expansion for the partition function 

n 

3f(n) = g-Mn|l U z ^ ( 3 - 2 °) 

X 1 i ■>■ ■■ > x n Jcomp 

where the sum is now over families of compatibles contours of the same class 
and 

z(r p ) = o(T p ) q ^0^^ (3.21) 
21 



with m 7^ p (and ||r p || is the number of lines of supp T p ). 

To control the behavior of our system, we need to show Peierls condition, 
that means that g(T p ) g e pH rP l! has good decaying properties with respect to 
the number of incorrect points |T P | of suppT p . In fact we shall consider 
the modified Peierls condition introduced in ref. pi] where e p is replaced by 
e = min (ef , ef Q ) . Let 

2 2d - 1 q ~WJ=T) + 24cu 3 d q~^~^ b A n r (3.22) 



1 - 6cu^q~^~ Pb ) 

2 

where c = 8e(e — l)(v / 3 + 2)e^+ 1 and is d = (2d) 2 . We have the following 

Proposition 3.3 (Peierls condition) Let Q C be a finite connected 
set of lines, assume that [3b < \, and 6cz/|g"d +/3t < 1 ; then for all j3 s G M : 



E 



e(T p) | gSFII < e -r||s|| (3 23) 



TP:suppYP=S 

where \\S\\ is the number of lines of S. 

The proof is postponed to the Appendix. 

The recurrence equations of Lemma EOl plav in our case the role of Lemma 
2.3 from Ref. (321 • Together with the Peierls estimates 13.31 they allow to 
study the states invariant under horizontal translation (HTIS) of the hydra 
system by applying Pirogov-Sinai theory. Here we do not need any extension 
[TUj [H] developed for interacting contours. Actually we have non interacting 
contours and only a standard form of the theory is required. We present it 
in the next subsection. We chose Zahradnik's approach [21] with the notion 
of truncated contour models including some improvements from [20J. We are 
dealing with complex activities, a situation analyzed in Refs. jHS] , jl] and we 
mainly follow the presentation given in that last reference. 

3.3 Diagram of horizontal translation invariant states 

To state our result, we first define the functional 

K P {S)= £ z(Tn (3.24) 

TP:supp FP=S 
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Consider the partition function E£ (Q) (|3.2(J|) and for a compatible family 
{T^, r^} comp of p-contours, denote by Si, S n their respective supports. 
By summing over all contours with the same support this partition function 
can be written as the partition function of a gas of polymers S with activity 
K P (S) = ^ ^(r p ) interacting through hard-core exclusion potential: 

rP:suppl>=S 



-dil 
J 7-' 



Here {Si, S n } comp denotes compatible families of polymers, that is 
d(Si, Sj) > 1 for every two Si and Sj in the family: recall that by definitions 
of contours a polymer S is a set of lines whose graph that joins all the points 
of the lines of S at distance d(i,j) < 1 is connected. 

Next we define auxiliary models, namely the truncated contour models 
associated to each boundary conditions. 

Definition 3.4 A truncated contour functional is defined as 

R , = (K P (S) if \\K P (S)\\ < e-^W 
p [ otherwise 

where \\K p (S)\\ = J2rp-su PP rp=s \ Z (^ P )\ > an do> > is some positive parameter 
to be chosen later (see Theorem \3.b\ below). 

Definition 3.5 The collection {S, p} of all p- contours T p with support 
suppT p = S is called stable if 

\\K P (S)\\ <e- a ^ (3.27) 

i.e. ifK p (S)=K' p (S). 

We define the truncated partition function H' (fi) as the partition function 
obtained from S^ 11 ^) by leaving out unstable collections of contours, namely 



n 

{r 1 ,... ) r^}comp 



Q 



-•pimii J2 U k p^) ( 3 - 28 ) 



{Si,...,S„} comp i—1 
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Here the sum goes over compatible families of stable collections of contours. 
Let 

^ = -Mmp|rln^(fi) (3.29) 

be the metastable free energy of the truncated partition function E' p (Q). 

For a large enough, the thermodynamic limit ()3.29j) can be controlled by 
a convergent cluster expansion. We conclude the existence of h p , together 
with the bounds 

e -«-|9.n| < sj^eMinil < e «-|8.n| (330) 
\h p — e p In q\ < ne~ a (3.31) 

shown below. Here k = k c \ (x') 2 where k c \ = v/ ^ +3 e^ g + 1 is the cluster con- 
stant [20] and re' = 3 d - 1 - 1; d s O = dtt nL in the way defined in Subsection 
PI 

Theorem 3.6 Assume /3b < l/d and q is large enough so that e~ a = 
e -r+2n +3 < then there exists a unique fit = + 0(e _r ) such that : 

(i) forP s = Pl 

s*(n) = s;(fi) 

/or &o£A boundary conditions p =f and p =of , and the free energy of 
the hydra model is given by gf+lnq = hf = hf . 

(ii) for(3 s <Pl 

s*(n) = s$(n) 

and gf + In q = hf < hf 

(iii) for(3 s >(3l 

and aj+ lng = /ij < /iy 

Let us introduce the Gibbs states ( • ) p , associated to the partition func- 
tions Sp 11 . Theorems above show also that at (3 S = /3* 

(n 6 ) f = 0(e^) 

fo _ i rk(„-r\ 



n b y° = l-0(e" 
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for any bond b of the boundary layer while 

(n b ,) { =0{e~T) , , 

for any bond b' between the boundary layer and the first layer. This is 
because, with free (f) boundary conditions, bonds b of the boundary layer 
and bonds b' of between the boundary and first layer are occupied only if 
there is a free contour surrounding them and that the correlation functions 
are controlled by the contour model cluster expansion. With free-ordered (fo) 
boundary conditions, bonds b are empty and bonds b' are occupied only if 
there is a fo-contour surrounding it and again the correlations are controlled 
by cluster expansion. 

This shows in particular that the derivative of the free energy <?f 

with respect to the surface coupling constant K is discontinuous near K = 

Obviously, the first relations of ()3.32j) and ()3.33j) hold true for any (3 S < /?* 
while the second relations of ()3.32|) and ()3.33j) hold true for any /3 S > 
Notice also that the estimate (|3.33|) on the mean values of the bonds b' 
between the boundary and first layer can be extended to the bonds of the 
bulk. 

Outline of proof of Theorem 13.61 

The first step is to prove the inequalities ()3.30|) and (|3.31|) . To exponen- 
tiate the partition function 

n 

Z p (Q) = q^E' p (Q) = IIW ( 3 - 34 ) 

{Si,...,5„} [ . omp i=l 



we define the truncated functional $ T associated to K' 

^ T ^) = Y^syUK(s) x(s) (3.35) 

defined on the multi-indexes X associated to the polymers: a multi-index is 
a function from the set of polymers into the set of non negative integers, and 
a(X) is defined as in (j2.11j) . 

The number of polymers S with number of lines ||S|| = n and containing 
a given line can be bounded by v n where v = (3 d_1 — l) 2 . Actually, There 
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exists at most 3 d_1 — 1 lines at distance 1 of a given line and for any connected 
set of lines there exists a graph going through every line at most twice. 
As a result of the standard cluster expansion (2II1I27|, we S e ^ f° r Ke ~ a < 1 

Z P (Q) = exp|^$ T (X)| (3.36) 

= ex P Uii y „^„-r iisu PP xnoii 

(3.37) 

Here, the sum in ()3.36|) is over multi-indexes whose support supp X = 
{S : X(S) > 1} belongs to Q, the first sum in ()3.37j) is over all multi-indexes 
X whose support contains a given line L, and the second sum in ()3.37j) 
runs over multi-indexes whose support intersects both Q and its complement 

L M \a 

The series Yl \& T (X) \ is absolutely convergent and satisfies the bound 

\® T (X)\<Ke~ a (3.38) 

X:snppX sL 

that immediately gives ()3.30|) and ()3.3H) . 

The second step closely follows the Borgs Imbrie paper 
We put 

a„ = h p — min h m (3.39) 

m 

The boundary condition p is then called stable if a p = 0. Our task is now 
to show that if the boundary condition p is stable, then all collections of 
p-contours are stable implying that Hp(f2) coincides with Ep ll (Q). 

We notice that when a p < 1 and the hypotheses of Proposition 13.31 hold, 
then the condition 

Sf (Int m P0 " l '- 4UJ 

for all contours T p with support S ensures that the collection {S, p} is stable, 
provided 

e~ a = e - T+2K ' +3 < - (3.41) 

K 

Indeed since 

z(Tn = g(r^)g^" rP l'g^-^l' rP l' 3IS! ntm rl (3-42) 
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we get, by Proposition 13.31 using (e p — e)\nq < a p + 2ne a < 3 (see ()3.31Jl ) 
and the inequality |<9Jnt m P>| < k' \\Y p \\, that \\K P {S)\\ < e --!l s ll e ( 3 + 2K ')ll s ll. 
We then show inductively on diamf2 = maxij e n n L d(i,j), that 



(i) if a m diamf2 < 1, and a p = 0, then 



' — 'm 



(O) 



< a m \\n\\+2Ke- a \d s n\ 



(ii) if a p = 0, then 



=dil 



(n) 



< e 3Ke- a |a s a| 



(3.43) 



(3.44) 



(iii) if a m diamf2 < 1, then 



1 — 'm 



< (l+5«;e- a )|9 s n| 



(3.45) 



The proof is analog to that of Theorem 3.1 in |4| using our previous 
estimates and thus we omit it. It gives that for ft < 1/d and q is large 
enough so that e ~ T+2K +3 = e~ a < then if a p = all collections of p- 
contours are stable. 

For our last step, it remains to show that there exists a unique /3* = 

+ 0(e~ T ) such that: 



af = and af = 
Of = and cif > 
af > and Of = 



for ft = # 
for ft < # 
for ft > (31 



(3.46) 
(3.47) 
(3.48) 

p^lnZ p (fi) is a 

Lipschitz function of ft uniformly in f2 with a small Lipschitz constant (when 
r is large, one has a good control on the one sided derivative ■g£^-K' p (S) of 
the activities associated to the partition functions Z p (Q)). 



This is a consequence of the fact that the free energy 



Lemma 3.7 Assume that ne a < 4 < ^r, £/ien 

ft' — k' 7 

1 9 



I^Hlng^ft 



lnZ p (ft) 



< (d 



a r K,e 



a r K,e 



(3.49) 



where a r 



3 9c/8 



h',7 
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Proof. 

By virtue of relation ([3.21)1 one gets 



d 

Ws 



K P (S) 



TP:supp TP=S 



d_ 5^(Int ro n 



(3.50) 



Since 



de p 



( ,r §g = -(d - 1) (see <KTT21t ). the bound (TQ5|) on the 



< 



first term inside the parenthesis gives that 

||r p || (d — 1) lng. This inequality together with the bound 1)3.14)1 on the third 
term inside the parenthesis leads to 







K P (S) 



< Yl [ll rP H {d-l)hiq+ ||Int m H (d-l)\nq\ \z(T p )\ 

TP;supp TP=S 

< [\\V(S)\\(d~l)\nq]e~ a ^ (3.51) 
where V(S) = V(T P ) =Int T p U P\ 

Taking into account the fact that either the one-sided derivative g£^K' p (S) 

d jy-i I Q\ ^ d TS I 0\ „, — +u„ l„f^- U„„ J „; J f /T5~l 



or gi^K'^S) < ^- s K p (S) we get that the left-hand side of (gUj) is bounded 



as 



1 5 



■lnZp(n) 



S-.SiL 



d 



d ± l3 i 



^ {Si,... 



{Si,...,S n } c 



<(d-l)\nqJ2 \\V(S)\\e' a ^expU ^ $ r (X) 

S:SaI< ^ X:suppXriS=0 

< (d- l)lng ^ 11^(5) || e -«llsil e (i+«')«e- Q lls|l ( 3 . 52 ) 



is over all families {Si, 5^}^ compatible with 5* 

and S = S U {L : S 1 ) = 1} is the set of lines at distance less or equal to 
1 from S. 

The last inequality leads to the lemma by inserting the bound 
\\V(S)\\ < ||S||diamS< \\S\\ 2 <e^ 
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and using the fact that the number of polymers with number of lines \\S\\ = n 
and containing a given line can be bounded by v n = (V) 2n . ■ 

To prove relations f)3.46H3.4H^ . we consider the energy difference ^Ej 1 = 
P s — -pp As a function of f3 s it is obviously increasing, negative for (3 a < -pr, 
positive for (3 S > ^rr, moreover it intersects the horizontal coordinate axis 
only at one point (3 S = ■j-r. 

Since ^"'"^ differs only from e ^'~ 1 ef by a Lipschitz function of j3 si namely 
the function lim^tL \\n\\(d-i)inq [l n -2-f(^) ~~ ^2f Q (fl)], it will satisfy the same 
properties (the intersecting point slightly changed) provided the Lipschitz 
constant of this last function is sufficiently small. Indeed, 

" r "~" 1 ' ' lini „ ' - [lnZ fo (fl)-lnZ f (fl)] (3.53) 



(d-l)lng d-l n?L ||f2|| (d-l)lng 

By virtue of Lemma 13.71 the needed condition on the Lipschitz constant is 
^ l-a^e-" < ^ since the Lipschitz constant of e ^J 1 ef is 1. This last inequality 
requires a c Ke~ a < 1/3 and is actually fulfilled under the hypotheses of the 
Theorem 13.61 by noticing that the geometric constant k' = 3 d-1 — 1 > 8, 
these hypotheses imply a c ne~ a < a c c/ ' k' < a c c/8 where the last term is less 
than 1/3 for c < 0.7. ■ 



Acknowledgments 



The authors thank S. Miracle-Sole, S. Shlosman, and V. Zagrebnov for helpful 
discussions. L.L. acknowledges the warm hospitality and financial support 
of the Centre de Physique Theorique. 



29 



Appendix: Proof of Proposition 13.3 



We begin the proof by considering contours T = {supp T, X r , F r } (where 
{X r , y r } is the configuration having T as unique contour) without polymers, 
i.e. Y r = 0. We let Xf = X r C)B(lL ) be the bonds of X r that belong to the 
boundary layer and X b = X T \ Xf . A site i G Lo of the boundary layer will 
be called regular if the bonds of its envelope that belong to the boundary 
layer are either all empty or all occupied. It will be called irregular otherwise 
and we will denote by /o(r) the set of irregular sites of the contour V. Namely 

i (r) = {% g l : i < \E(i) nxf| < 2(d - 1) - 1}. 

Lemma A.l 

l^)|^imi< g -^f-(^)K p l (A.i) 

Proof. By Lemma f3. 21 and definition ()3.9|) . one has 

g (T p )q^ = q -( H su PP M x )-e\\r p W) 

= g £ iesupp rPn S (L ) ^\E(i)nX s \+MX b \~S(X) + C(X)- Xp +e\\rP\\ / A 2 x 



Here, Xf = 0, Xo — 1) an d to simplify formulae we put hereafter X, X s , and 
X b instead of X rP ,X^ P , and X[ P . By using the relations 

S(X) = S(X S ) + S(X b U X s ) - S(X S ) (A.3) 
C(X) = C(X S ) + C(X b U X s ) - C{X S ) (A.4) 



we get 
where 



0(rP)gSP*ll < q As(VP) + A b (T*) (A 5) 



A s (T p )= (j\E(i)nX s \-\inS(X)\+zj+C(X s ) (A.6) 

A 6 (r p ) = (3 b \X b \ - S(X b U X.) + S(X S ) + C(X b U X.) - C(X S ) (A.7) 

Let us first bound A s (r p ). 

Consider an irregular site (of the boundary layer). Then by definition: 
1 < \E(i)nX s \ < 2(d - 1) - 1. Since e= min{e f = 0, e fo = 1 - (d - 1)/3 S }, 
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we have, 



y \E{i) n X s \ - \i n + e = y n X s | - 1 + e 

S 2^1)l £ W nX -|- 1+£ 
nl s | -2(d- i; 



1-e 



2(d- 1) 



|£(»)nx,|-2(d-i) 
" 2(d=T) (A ' 8) 

where we used successively /3 S < (1 — e)/(d — 1) and e < 

Next, we take into account that the number of connected components 
may be bounded as (c.f. e.g. [T§1 122J): 

i:l<\E(i)nXs\<d-l 

When 1 < \E(i) n X s | < 2(d - 1) - 1, we infer 

\E(i) nx s \-2(d- i) x(i < I^W nx,|<d-i; 



< — — - (A.10) 



2(d-l) 2l s « n ^l ~2(d-l 

so that each irregular site i of the boundary layer provides at most a contri- 
bution 2 (d-i) ^° ^ ne right-hand side of ()A.6|) giving 

g A s (rf) < 5 ^317 (A. 11) 

Consider now, the quantity A^(r p ). One has 

q A b (rp) < Y^ q i 3 b\s\-s(sux s )+s(x B )+c(Sux s )-c(x s ) (A. 12) 

where the product runs over connected component of X^. By the same ge- 
ometrical observation as in the proof of Theorem 12.11 one easily gets the 
inequality 

^--S{8U X.) + S{X S ) + C{5 U X a ) - C{X S ) < (A. 13) 
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that gives 

q A b (rp) < q -{\-M\x b \ (A. 14) 

which in turn implies the lemma by taking into account (|A.5|) and (|A.11|) . 
Notice that by ()A.2|) one has the following bound 

° \yiq(T 1} ) < \\Y P \\ (d- l)lng (A.15) 



Considering still contours T = |supr,X r ,y r } without polymers, i.e. 
(Y r = 0) we have the 



Lemma A. 2 Assume that f3 b < \, and q<i & b > 2v d = 2{2d) 2 , th 



en 



qF>) ^ l|rP)l < (2 2d - l q-^ + Su d q-(^)) m l — — 

TP:suppTP=S 1 — 2f d g V<* h ) 

(A.16) 

which shows that, whenever q is large enough, the Peierls condition holds 
true for the class of contours without polymers. 

Proof. First, observe that for contours T p with support suppT p = S 
and number of irregular sites of the boundary layer |/ (r p )| = k one has 
\Xb\ = + ••■ + \5 m \ > \\S\\ _ k. Therefore, 

rP:sup P rP=5 0<fe<||5|| rP:\I (TP)\=k 



< Y ^ll' 5 l| N J 2 2(d-l)fe 2 ||S||-fc g -^ TJ 

L.S'W Oil V ^ / 



0<fc<||5|| 



X 



lib 

e e e n*-^ 1 * 1 



ra<||S|| 8i3ii,...,8 n 3i n j=l 

icSSnLo |5i|+...+[5n[>|[-Sf||-A! 



(A.17) 



Here the binomial coefficient ( f ) bounds the choice of irregular sites of 
the boundary layer while the factor 2 2 ( d-1 ) fc 2" 5 "~ fc bounds the numbers of 
contours with 1 15*11 lines and k irregular sites. Then 
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E e(r*) 5 ""< E ("f") 

rp : sup P rp=5 o<fc<||5j| ^ 7 



22(d-l)fc2l|g||-fc - a^,!) 



x 



|s „s 



e "i" e n k H) < a ' i8 > 



n<||5|| ^ ' mi + ...+m n >||S||-fc J=l 



Here the binomial coefficient ("f") bounds the choice for the components 
Si, ...,S n of X b to hit the boundary layer at ii, ...,i n . The above inequality 
yields 

,™n.-.c 0<fe<||S|| ^ ' 

e (»f ) e (*«-«-*>: 

<\\S\\ v 7 m>||S||-fc 

< ('"' S '" N )2 2 ^- 1 ) fc 2ll s ll- fe g-^T 

n^k^ll oil V / 



rp : sup P rp=5 o<fc<||s| 

X 

n<||S|| V 7 m>\\S\\-k 



0<fe<||5!| 



x (2^-(3-^)) (A.19) 

that gives the inequality of the lemma. ■ 

We now turn to the general case of contours with non empty polymers 
and first give a bound on the activity ip (7) of polymers. 

Lemma A. 3 Assume that /3& < | ; and cv 2 d q^^^ b < 1 c = 8e(e — 
lXv 7 ^ + 2)e^fe and i/ d = (2d) 2 , then 

IV (7)1 < (a^?-^ +A ) H (A.20) 

Proof. Let us first recall the definition (|3.6|) : ^(7) = J2a-su PP a=j 
where the weights of aggregates are defined by (see (|3.3|) and ()3.4j) j: uj(A) = 

n 7eA e"* (7) -l. By Theorem Owe know that |$( 7 )| < (ecQP d q-^ +Pb \ H (< 



1) for q large enough. Since for any |x| < 1, |e x — 1| < (e — 1) we have 
1^(7)1 = |e~* {7) - 1| < (e - 1) |$( 7 )| < ((e - l)ec u d q-^ 



-<r|7| 



e 

(A.21) 
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Then, 



E w)i=E E ni»(-» 

.4:suppyl=7 n>l 7l>— >7n j=l 

supp{7i,...,7„}=7 



n 



<E 2W e ri e ^ i7ji 

n>l 71 =*i,..., 7„3&„ i=l 

supp{7i,...,7„}=7 

n 

<E 2hi e n^-n 

n>l mi,...,m n J= i 

mi+...+m n >|7| 

n 

^E E II( 2 ^~ CT n ( A - 22 ) 

n>l mi,...,m n j=l 
m 1 + ...+m„>\-y\ 

Here, we used as in the proof of Theorem 12 . II that the number of polymers of 
length m containing a given bond or a given vertex is less than z/^ m ; the term 
2^ bounds the combinatoric choice of the bonds or vertices bj e jj, because 
7 being connected, it contains n — 1 such intersecting cells (see 

We put k = rrii + ••• + rn n and notice that there are at most ( _J such 
numbers to get 

E k»i= E E(„*iW~)* 

A:suppvl=7 l<n</cfc>|7| ^ ' 



k>\j\ k>\-y\ 

< ( -^<T* + M hl (A.23) 

provided that §i / Jq , ~^ +/3b < 1. The lemma then follows by assuming that 

fi&-* +A < | ■ 

We finally turn to the 
Proof of Proposition 13.31 

Consider a contour T = {supp T, X r , E r } and as above the decomposition 
X r = Xf U X[ where Xf = X r n S(L ) is the set of bonds of X r on the 
boundary layer and X^ = X T \ Xg. Consider also the union = X^ U Y T . 
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Notice that the set Z = U Zj, is a family of hydras and there are at most 
3\ z b\ contours corresponding to this family: this is because a bond in Z b may 
be occupied either by X^ or by Y r or by both. Let 

\g(Z)\= l^ F )l (A.24) 

r-.xluY r =z 

The above remark on the number of contours associated to Z and Lemma 
,~3l implies 



|g(p) \q &m <q ^3sup{g~(3~ A ) >C i^ 9 -3+A| 



|z»l 



IMQI / 9 1 , R \ \ z b\ 
<g « 3ci/ d V 3 A (A.25) 



The rest of the proof is then analog to that of Lemma IA.2I starting from 
Lemma lA. II and replacing VdQ~ 1+l3b by 3cu^q~d +/3b . It gives 



l^( r? ) I ^"" rP " - (2 2d ~V^ + 24cz/Jg 



\\s\\ l 



(A.26) 

provided 6cu^q~d + ^ b < 1 and ends the proof of the proposition. ■ 
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Figure captions 



1. The diagram of ground states. 

2. A hydra, in two dimensions (a dimension not considered in this paper), 
with 2 feet (components of full lines), 5 necks (dotted lines), and 3 
heads (components of dashed lines). 



, 1 _ _ _ J . 



